Excitations of single-beauty hadrons 
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In this work we study the predominantly orbital and radial excitations of hadrons containing a 
single heavy quark. We present meson and baryon mass splittings and ratios of meson decay con- 
stants (e.g., fBs/fs and fgi / fs^) resulting from quenched and dynamical two-flavor configurations. 
Light quarks are simulated using the chirally improved (CI) lattice Dirac operator at valence masses 
as light as ~ 350 MeV. The heavy quark is approximated by a static propagator, appropriate 
for the b quark on our lattices (l/a~l — 2 GeV). We also include some preliminary calculations of 
the Oil/mq) kinetic corrections to the states, showing, in the process, a viable way of applying the 
variational method to three-point functions involving excited states. We compare our results with 
recent experimental findings. 

PACS numbers: 11.15.Ha, 12.38.Gc, 12.39.Hg, 14.20.Mr, 14.40.Nd 



I. INTRODUCTION 

A. Singly beautiful hadrons: Recent progress 

With one heavy valence quark and otherwise light de- 
grees of freedom, single beauty (bottom, 6) mesons and 
baryons are the QCD analogs of hydrogen and helium 
atoms. In the past few years, there has been much 
progress in determining the properties of excited states 
of such hadrons. The DO and CDF experiments at Fer- 
milab have recently observed the "orbitally excited" Bi , 
Bl j,] and Bsi, i?*2 mesons as well as the S^, Y.*^ 



and rifj baryons [5;]. 



b ' 

There have also been 



advances in lattice calculations of the mass 



ectrum of 

excited B mesons ^ 

and h baryons 0, [3113 l2Cll. l2l|. as well as calculations 
involving hadronic decay of the Bq [2^ . 

There are certainly daunting problems involved in 
the experimental efforts, and likewise, extracting excited 
hadrons is a real challenge on the lattice. Excited states 
appear as subleading, decaying exponentials in noisy 
two-point correlators and separating a given state from 
lower- and higher-lying states becomes a highly non- 
trivial problem since the former dominate at large Eu- 
clidean times, while the latter tend to be close in mass 
and therefore give similar contributions over the same re- 
gion. There are a number of methods which have been 
devised to handle exactly this problem. In this paper, 
we focus on just one of these, the variational method 
[ll,!!!] (see also [H,!!!,!!!!!!): From matrices of two- 
point correlators we determine meson and baryon mass 
splittings and meson decay constants; we also expand 
the method to handle three-point functions, in order to 
calculate 0{\/mb) kinetic corrections to our static ap- 
proximation of the b quark. Preliminary mass s plit tings 
from this ongoing study were reported in Refs. p^. 
more recent results appeared in [29|. 



This paper is laid out as follows: In Sec.|TT]we discuss 
the variational method and how we use it to extract the 
masses, couplings, and kinetic corrections of our heavy- 
light(-light) hadrons. The details behind our lattice cal- 
culation, including our basis of operators, are presented 
in Sec. IIIII In Sec. IIVI we present our analysis of each of 
the items mentioned in HIl and in Sec.|V]we discuss our 
results for each of the mesons and baryons. We conclude 
and present an outlook in Sec. IVII 



B. Operators versus states: Not just semantics 



Before moving on, we believe it necessary to point out 
that the words "orbital," "radial," and "gluonic excita- 
tions" should be viewed as only convenient terminology 
since each describes only one (albeit perhaps the most 
likely) of the configurations of the constituents in a par- 
ticular excited state. In QCD, mixtures of many such 
configurations are present in any hadronic state [s^, [3l| 
and their clear separation only begins to make sense in 
the heavy-quark limit (3^.[33|. The same applies to many 
similar terms, like "hybrid" and "P-wave," and to any 
definition which tries to fix the number of constituents 
since higher Fock components also represent ever-present 
admixtures. What this means is that one must take care 
in distinguishing the properties of operators (or interpo- 
lators) from those of states. Ultimately, it is the covariant 
nature of the derivative (the momentum operator) which 
ensures that orbital, radial, and gluonic excitations are 
not uniquely defined. 
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II. METHODOLOGY 



A. Masses 



Since we are using a variational method, we need to 
choose a large enough basis of operators to create and 
annihilate the states of interest. We discuss the details 
of our choice of basis in the next section, but assuming 
for the moment that we have these, the next step is to 
create all cross-correlators between such operators: 



C(% = (0|O.(i) O](0)|0) 



(1) 



We may write each of these correlators as a spectral sum 
of energy eigenstates 



cm 



(n) 

where the represents the overlap of the ith operator 
with the nth state 



E(n) (n)* —tE^"^^ 

n=l 



(2) 



,(") 



:o\o,\H. 



We use the notation IH^"^) for the nth hadronic state 
containing a single static quark (Q). 

Solving the generalized eigenvalue problem 



Q . 



(3) 



then leads to eigenvalues of the form [23| 



(4) 



X{t,to) 



(a) 



l + 0{e 



-(t-*o)A(° 



(5) 

where it is the sum of all terms which equals 1 at t — to 
and the quantity A^"^ represents the minimum energy 
difference to an adjacent state (to the first state beyond 
the basis dimension for to < t < 2to [ll]). So with a 
clever choice of basis (which can reduce the amplitudes 
of the correction terms) and for large enough values of 
t (ideally, t ^ to), each eigenvalue yields the mass of a 
single state. The eigenvectors ■(/'^"^ can be used to provide 
information about the couplings to the state a (see Ref. 
[26t and the next subsection). We do not indicate an 
explicit ^-dependence of the eigenvectors since, ideally, 
there would not be any; however, one should realize that, 
due to higher-order corrections (see Eq. pH]) below) and 
statistical fluctuations, there is in fact one. 

Following heavy quark effective theory (HQET; see, 
e.g., Ref. |3J|), the masses of hadrons containing a single 
heavy quark may be written in terms of an expansion in 
I/tuq. Accordingly, the heavy- light masses are 

= mQ+ -Eo + 0(fc(")V™Q) , (6) 

where A^"^ is the static-light binding energy of the nth 
state and fc'-"-' denotes the residual momentum of the 



heavy quark in the nth state [3^ (for the ground state, 
fc(i) Aqcd)- Note the shift Eo, indicating that we do 
not have a priori a physically meaningful, zero-point en- 
ergy of the static quark. In fact, appropriate definitions 
of mg and Eq must be chosen to avoid power divergences 
and renormalon ambiguities [36| in any determination of 
nij or A(i). We may, however, use the lattice energy 
differences to determine the mass splittings 

^("') _ ^(") = a'"'-* — A*-"' 

= Af("') - m(") -h 

-fO(A:("')VmQ,A;(")VmQ), (7) 

including those between hadrons of different flavor, e.g., 
for Si*^ [13, 

£;(!) _£;(!) - A(i)-A(^^ 

= Mgi,) -Mbm + 

+ 0{ki^^ym,,k^^-)ym,). (8) 

In particular, we use this quantity to set the physical 
strange quark mass (see Sec. IIIIE)t . 



B. Couplings 

If we work at large enough values of t, then we may 
trust that, up to higher-order corrections, each side of 
Eq. ^ describes a single mass eigenstate. Working with 
the left side of Eq. (g]), we find 



(9) 



The solution of the generalized eigenvalue problem then 
ensures that 

J2 v^;^*^f^ = <5("'") + 0(e-(*-*°)^*°' ) , (10) 

where A*^"^ is an arbitrary normalization which one can 
easily remove (see below). Therefore, for large t values, 
we can write 12611 



(11) 



We can also use two eigenvectors, giving 



i=l j=l 

We may then construct and fit the following ratio to 
determine the couplings: 



R{t) 



(a) 



(q) (q)* _t_E(°) 



(13) 



3 



Ratios of different couplings to the same mass eigenstate 
[2al are even easier: 



,(") 

J. 

,(") 



(a) 



(14) 



Since we are working with the static approximation 
for the heavy quark, the local vector (O^ = qjiQ) opera- 
tor overlaps are related to the pseudoscalar (PS) decay 
constants via [3J| 



In order to cancel renormalization factors and coefficients 
matching HQET to QCD, we deal with ratios of over- 
laps. So, for example, the ratio of meson decay constants 
fs' I rnay be extracted from the ruq — mg point of 



f(2) 

J PS 

J PS 



„W Y M^(., + - Ed)) 

+ 0(fc(2)Vm,,fcW2/^^) ^ (16) 



where = 5400 MeV and the and E^"^ come 



from fits to Eq. (|13p for a — 1 and 2. Similarly, we can 
find fsJfB from 



J PS 



,(1) 



I Mr,,, 



B(') 



J PS 



„(i) 



m,„=m^d 



irLa—m-,, 



where M^co = 5314 MeV. 



(17) 



C. Kinetic corrections 

In order to go beyond the static approximation for 
the heavy quark, we begin calculations of the O^I/uiq) 
kinetic corrections to the static-light states. The new 
term in the Lagrangian is just the kinetic energy of 
the heavy quark and the corresponding operator is the 
lattice-discretized Laplacian. The relevant matrix ele- 
ments are then 



.(n,n') 



(«') • 



(18) 



What we call e^.i) corresponds to the negative of the 
usual HQET parameter Ai. From the e("'" ) values, the 
energy shifts {n — n') and static- light state mixings [n ^ 
n') can be calculated: 



2m, 



-Z{e 



(n,n) 



eo) , 



(19) 



Sv, 



(") 



1 



2m 



E 



Z{e 



{n' ,n) 



^o) (n 



(20) 



The subtraction of eg represents the need to remove the 
ultraviolet (power-law) divergence in the lattice version 
of the QD^Q matrix element and the Z denotes the cor- 
responding renormalization factor for this operator [36j . 
For a more complete, future study we will need to calcu- 
late these quantities; however, our focus here is the tech- 
nique we use to find the bare corrections for n, n' > 1. So 
for now we report differences of the e("'" ) to cancel the 
first quantity and use a tadpole-improved js^ version of 
(dividing the spatial links used in the operator by uq) 
to try to tame the second. 

To extract the bare corrections, the quantities of inter- 
est are the lattice three-point functions 

T{t,t% = -{0\O,{t) QD^Qit') O]{0)\0) . (21) 

The trouble with applying the variational method di- 
rectly to this new matrix of correlators is that all entries 
involve a valence quark participating in the current in- 
sertion QD^Q. This interaction therefore excites its own 
tower of states in the correlators Tij , much like the source 
(Oj) and sink operators (Oi) do. 

One possible solution (yet an approximate one; see dis- 
cussion below) is to view the three-point functions as 
taking part in two variational problems: 

C{t-t% = {0\O,{t)O}{t')\0) 

C{t% = (0|O.(t')O|(0)|0>. (22) 

Solving the two corresponding generalized eigenvalue 
equations, 

c(i-i')V^^"^ = A(t-t',t^,-t')^"^c'(t;,-i')V^^"^ 

Cit')$^^^ = X{t',t„)^^^C{to)$^^^ , (23) 

then gives sets of eigenvectors which one can use to 
project the states of interest on either side of the inter- 
action: 

^^#)*r(M').0f 

z=l j=l 

_^(a)*g(a,/3)^(,3)g-(t-t')S(°'g-t'-B<« _ (24) 

With a large ensemble and all-to-all quark propagators, 
translational invariance should strictly apply and one can 
use just one set of eigenvectors V''"' — <t'^°'\ referenced 
by just the time difference At = t' ot t — t' (this we do.) 
One may then form ratios to cancel the exponentials and 
some of the overlap factors: 



n 



kl 



£(",/3) 

(a)* (/3) ■ 



{13) 



(25) 
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Similarly, it is possible to contract the open k, I indices 
with the corresponding eigenvector components and end 
up with A'^°'^A^^''* in the denominator. Combining this 
with earlier ratios for determining the couplings, we ob- 
tain (assuming real overlaps) 

<-'^Vi?(t-t')^^i?(^')f 

_ g(a,/3) g-[(t-t')-E<"'+t'B<^']/2 _ (26) 

We can fit such a quantity directly for the bare kinetic 
corrections e^"'''^. (This technique and first results were 
presented in [2§|-) 

This shows how the variational method may be used 
to extract matrix elements and form factors (including 
transitions) involving excited hadronic states. The fact 
that we are dealing with ©(I/toq) corrections to HQET 
states is irrelevant. The procedure can be easily general- 
ized to, say, an electroweak current applied to full QCD 
states. 

The salient feature of our analysis is that the three- 
point functions are treated in the context of two varia- 
tional problems: one before the interaction and one after. 
It is conceivable, however, that a chosen basis of source 
and sink operators may not be robust enough to well ap- 
proximate the intrusion of the probing current. In other 
words, the current itself must be viewed as a linearly in- 
dependent operator exposing the higher dimensionality 
of the problem at hand (with lower lying states suddenly 
"reappearing" in a given eigenvalue channel a, perhaps 
first in the form of seemingly additional noise). Indeed, 
in the limit of large statistics, this should be seen with 
any finite basis in which the interaction is not included. 
To be completely accurate, the current, together with 
the evolution to the sink or from the source, has to be 
added to the basis as an additional operator. Such com- 
putations may be prohibitively expensive since the new 
diagonal elements of such an expanded correlator matrix 
would be four-point functions. 

One may also try a more simplified version of Eq. (j25p 
for the case oi a = /3: 

E,EA"'*c{t),4'^^ 

« , (27) 

where tp^"^ is now the eigenvector from the solution of the 
generalized eigenvalue problem for C{t), Eq. (jH). As long 
as t' and t — t' are large enough, using these eigenvectors 
or the ones from Eqs. (|^ should give consistent results. 

III. LATTICE CALCULATION DETAILS 

A. Configuration details 

We perform runs on a combination of both quenched 
{Nf = 0) and dynamical {Nf = 2) ensembles The 



gauge action is one-loop improved Liischer-Weisz ^Ad\ and 
we use (for both sea and valence quarks) chirally im- 
proved (CI) fermions [4l| . which have been shown to have 
good chiral properties We set the scale of the lat- 
tice spacing via the static-quark potential ji^ , using the 
value ro — 0.49(1) fm Except where noted other- 

wise, the links of the quenched configurations are Hyp 
smeared (45j and the dynamical ones are Stout smeared 
p^ . The details of the configurations can be found in 
Table H 

In Table |TT] we outline the basic results of a sim- 
ple test of our choice of a physical value for tq. Us- 
ing tadpole-improved, 0(?;^)-accurate lattice NRQCD 
(with unsmeared gauge links; see Ref. for details of 
the action) we create heavy quarkonium correlators at 
zero and finite momentum. We use the spin-averaged 
IP — T splitting (440 MeV) to set the lattice spacing 
for these runs and the kinetic T mass to set the physical 
b quark mass. The results on the P — 4.65 dynamical 
lattice give 1/a = 1231(7) MeV, which, together with 
the value of rg/a from the static-quark potential, gives 
To — 0.5027(29) fm. This value is quite precise, but incor- 
porates partial quenching effects since the sea quarks are 
still quite heavy (M^^sea ~ 460 MeV). However, we are 
encouraged by its consistency with our "working value" 
of 0.49(1) fm, which we have chosen to be slightly lower 
than our previous value (0.50 fm in Refs. [Ij, [lg|) due 
to a seeming average trend in the results of the larger 
lattice community (see, e.g., Ref. |47j ) . 

B. Light-quark propagators 

In this subsection we briefly cover the improved es- 
timates of half-to-half propagators which we use for 
the light quarks. Details of this so-called domain- 
decomposition improvement (DDI) can be found in Ref. 

. The technique is inspired by earlier work for pseud- 
ofermion estimators Q, but represents a more tractable 
method for improved lattice fermions. 

If we decompose the lattice into two distinct domains, 
the lattice Dirac operator {A4) may be written in terms 
of submatrices 



TABLE I: Relevant lattice parameters. The scales are set 
using ro = 0.49(1) fm [ij. 



13 


N'^ X Nt 


Nconf 


a-' (MeV) M^.sea (MeV) L (fm) 


8.15 


20-^ 


X 40 


100 


1682(34) 


oo 2.34 


7.90 


16^ 


X 32 


100 


1360(28) 


00 2.32 


7.57 


12^ 


X 24 


200 


1010(21) 


oo 2.34 


4.65 


16^ 


X 32 


100 


1263(26) 


461(6)(9) 2.50 


5.20 


12^ 


X 24 


74 


1750(36) 


525(22)(11) 1.35 
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where Adu and A^22 connect sites within a region and 
M.12 and A^2i connect sites from the different regions. 
We may also write the propagator (V) in this form: 

The propagator between regions 1 and 2 can then be 
estimated using N random sources (x% « = 1, .., N) [l3l |: 

i 

^ E (-^n >'i2X^) (75T'75X2)^ = 

-^E*1'^'^ (30) 

The 2N solution vectors (^''j^^) are calculated and saved 
to later reconstruct the estimated light-quark propagator 
from anywhere in region 1 to anywhere in region 2. In 
this form we see that no sources are needed in region 1 
and those in region 2 should reach region 1 with one ap- 
plication of M . Since M is usually a sparse matrix (even 
for the CI operator), this greatly reduces the number of 
lattice sites which the random sources need to cover and 
therefore greatly improves the signal-to-noise ratio for 
long-distance correlators. 

We stick with the implementation as originally de- 
tailed in p^ . In particular, each domain contains half 
of the timeslices, the sources appear in both regions at 
the boundary edges (in order to obtain two estimates 
of 7^12), and we divide our random sources (originally, 
N — 12) into a factor four greater in number by making 
the spin summation explicit (spin "dilution" ) . 



TABLE II: Scale setting on the /3 = 4.65 (ro/a = 
3.1414(21)) dynamical lattice using 0{v'^) NRQCD bb corre- 
lators, AS(1P- T) = 440 MeV (to set a~^), and Mr = 9460 
MeV (to set arrib). Result: ro = 0.5027(29) fm, consistent 
with our more conservative use of ro = 0.49(1) fm through- 
out this work. 



anib 


state t/a aE 


a-' 


^is.kin 


3.2 


IS 20-25 0.44338(29) 








IP, spin- 








averaged 3-8 0.8004(21) 


1232(7) 






IS @ 


MeV 






p = 2TT/aNs 20-25 0.45369(30) 




9214 MeV 


3.6 


IS 20-25 0.42580(26) 








IP, spin- 








averaged 3-8 0.7851(21) 


1225(7) 






15 @ 








p = 2n/aNs 20-25 0.43486(27) 




10426 



3.28 1231(7) 9460 



C. Operator basis and correlator construction 



The local versions of the operators we use can be found 
in Table IIIII This is just the start of our basis and we 
expand it further by gauge-covariantly smearing [48l | the 
light-quark wavefunctions. 

In choosing how to smear the light-quark sources (and 
sinks), we work with a different philosophy than we em- 
ployed previously (see, e.g., Ref. [4§|). Here, rather than 
designing the sources, we simply make choices such that 
the resulting operators are drastically different from each 
other (having, e.g., relatively small overlap in the free 
case). As long as we do not make the sources too large 
or too small, this should ensure that the states of inter- 
est (the lower-lying ones) have various amounts of overlap 
with the different operators, hopefully making the states 
easier to resolve. 

Accordingly, for each set of meson quantum numbers, 
in addition to the local operator, we construct three 
smeared versions: one left as is, one with a subsequent 
gauge-covariant Laplacian applied, and the last with two 
such applications. For the baryons, we have two different 
local operators and we use smeared versions as well (so 
far, only smeared on the /3 = 8.15 lattice). This allows 
us to use up to a 4 X 4 basis for each desired (re- 
duced to 3 X 3 for the mesons on the dynamical lattices 
due to lower statistics). The parameters chosen for the 
smearings are shown in Table ITVl 

Given the operators of the basis, the correlator matrix 
can then be constructed. For example, for the static-light 



TABLE III: Heavy hadron operators. Shown are the total 
light-quark (plus glue) angular momentum j (lowest value 
in the continuum), the possible quantum numbers when 
including the heavy-quark spin, and the spin-averaged heavy- 
chiral (-strange) states [IJ- The last column displays the local 
version of the operator; color (a, b, c) and Lorentz indices are 
explicit, while spatial and spin indices are implicit. 



label j 




states 


operator 




(0,1)- 


(Bi*') 


Qa 75, i Qa 


P- k 


(0,1)+ 


-Bo.i (-Bso.i) 


Qa 7i(A)a6 


P^ 1 


(1,2) + 


r(-) |'r(*) ^ 
-"1(2) \^al(2)l 


Qa {jlDl - 72-D2)a6<76 


D± (f,f) (1,2,3)- 


„(*',*) try(*',*)\ 
^2(3) y^s2(3) ) 


Qa 'y^iDl- DDabQb 


Aq 


1 + 
2 


A6 (-) 


<^abc Qa Qb C75 qc 








if^abc Qa qb C7475 qc 




(1 1) + 
V 2 ' 2 / 




(■abc Qa qb Cji qc 



i^abc Qa qb Cj4^i qc 



6 



1 8 3 



r ' ' ' r ' ' ' r ' ' ' I 



1 8 



r ' ' ' r ' ' ' r ' ' ' I 



1 2 



FIG. 1: The setting of the strange quark mass on the /3 — 7.90 and (3 = 8.15 quenched lattices and the (3 = 4.65 dynamical 
lattice (from left to right, respectively). Plotted are the 15 lattice energies ro-E ' versus {r^M^f. The si'' - E^^l difference 
is set to the static-strange-static-light mass difference of 76.9 MeV, which is the 1/M^(,) — > linear extrapolation of the 
experimental values Af^c*) — — 86.8 MeV and M^(,) — Mj^f,) — 103.5 MeV. The burst at the left shows the linear chiral 

extrapolation (-E^^') and the two bursts at the right show the linear strange-mass interpolations (-EI^') using ro = 0.49 and 
0.50 fm. 



mesons we have 

= (0|O,(t) Ol(0)|0 



,(31) 



where the product of the (smeared) links in the time 
direction corresponds to the evolution of the (spinless) 
static quark, the lower-case operators Oj correspond to 
what remains of the Oj after integration over the fermion 
fields (Q, Q, q, q), and V^^ti^ is the estimated light- 
quark propagator (discussed in the previous subsection). 
So the lattice points x and x + tA must be within different 
domains and the ensemble average implied by the angled 
brackets also includes a division by the number of correla- 
tors of extent t which cross the boundary: min{t, T — t). 
For the static-light-light baryons, there is a similar ex- 
pression, but without the (l-|-74)/2 projection and where 



TABLE IV: Parameters for the quark source smearings. 
Listed are the number of smearing steps Nsm,i, the hopping 
weight factor k, and the number of Laplacians applied k for 
each set of configurations. Except where noted otherwise, the 
link smearing is also in place for the propagation of both the 
CI light quark and the static quark. 



f3 


link smearing 




8.15 


Hyp [45] 


(0,0,1,2) (0,16,24,32) 0.2 


7.90 


Hyp 


(0,0,1,2) (0,12,18,24) 0.2 


7.57 


Hyp 


(0,0,1,2) (0,8,12,16) 0.2 


4.65 


Stout [46| 


(0,0,1,2) (0,12,18,24) 0.2 


5.20 


Stout 


(0,0,1,2) (0,12,18,24) 0.2 



the (now, two) estimated light quarks are constructed 
without using the same random sources in the corre- 
lator contributions (i.e., V^yV'^y oc ^^^y ^i^^'^ 
has the restriction that i ^ j)- 

Before proceeding to solve the generalized eigenvalue 
problem, we check that the correlator matrices are real 
and symmetric (within errors). We then explicitly sym- 
metrize them and use a Cholesky decomposition of C{to) 
to solve Eqs. (U) and (p3)) . 



D. 



First look at the ground state: Physical strange 
quark mass 



To determine the physical strange quark mass, we 
consider spin-averaged experimental heavy-light meson 
masses. Up to ©(l/mg) in HQET, the average mass of 
the heavy-light ground-state spin multiplet is given by 
(using our notation) 



^(1) 

9 



£o) 



2m, 



Q 



(32) 



Since we want to know the strange quark mass, we then 
consider the following difference 



M 



HM 



(1,1) _ (l,l)^ 



2mr 



(33) 



To minimize the effect of the l/mq term on the right 
hand side, we linearly extrapolate the experimental num- 
bers [50| M^(.) - M^i,) = 103.5 MeV at l/M^,,., and 
M^(.) - Mbm = 86.8 MeV at l/Mgc,) to 1/Mhm -> 0. 
This gives M^(.) - MHi.) =76.9 MeV and we use this 
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number to set the value of (roMTr)^ corresponding to the 
physical strange quark mass. 

Figure [1] shows the ground-state static-light 5-wave 
masses versus the corresponding light-quark masses 
(oc M^) on three of the ensembles we use (see the next 
section for details surrounding the fits). The data points 
are consistent with a solely linear dependence over the re- 
gion of simulated quark masses and we therefore perform 
only linear extrapolations (to rriud) and interpolations (to 
iris). 

Given the quality of these results and the completely 
linear behavior they display, one can imagine using them 
as input in a determination of 7n7(2 GeV). We plan to 
pursue this in the near future. 

IV. ANALYSIS 

In this section we present our analysis of the static- 
light(-light) lattice energies and couplings. These lead to 
the corresponding results for the meson and baryon mass 
splittings and ratios of meson decay constants. Lastly, 
we show our first results for the bare kinetic corrections 
g(n,n ) ^Yie nS and 1P+ meson states on a subset of the 
quenched a « 0.15 fm (/3 — 7.90) configurations. 



the 25* states on the /3 = 4.65 dynamical lattice, where we 
use to/a = 1, but there we fit A*-^^ with a two-mass ansatz, 
trusting only the lower mass. Similar two-mass fits are 
also used for the 25" states on the smaller, finer, dynami- 
cal lattice, as well as for all of the ground-state baryons. 
For the remainder of the results, the bases seem to work 
well, suppressing the higher-order corrections (see Eqs. 
and pop ), and we can perform many single-mass fits. 
As mentioned previously (Sec. IIII Cp . most of the re- 
sults come from the full 4x4 basis. The exceptions 
are: the mesons on the dynamical lattices, where we find 
it beneficial to prune the basis to the first 3 operators 
(with the full basis, the solution of the eigenvalue prob- 
lem becomes unstable since, within the statistically lim- 
ited ensemble, the fourth operator adds more noise, but 
seemingly little new information); and the baryons on 
the P — 8.15 quenched lattice, where we only have re- 
sults from the smeared light-quark sources thus far (2x2 
basis) . 

In Tables IVl - IIXI we show details of our best fits for the 
IS", 2S, IP-, 1P+, 1D±, Aq, and S^*^ states at ah values 
of the simulated valence quark masses on each ensemble. 
(In the interest of some brevity, we do not present such 
detail for all states; there are, however, some 35, 2P, 
2D, and Ag results which appear in the effective-mass 
and chiral-extrapolation plots and later tables.) 



A. Masses 

1. Diagonalization, effective masses, and fits 

Constructing the correlator matrices and solving (or 
diagonalizing) the corresponding generalized eigenvalue 
problem as outlined in Sec. IIII C| we then have a set of 
eigenvalues X{t, io)^"^ and eigenvectors i/;^"^ (up to n = 4 
for the 4x4 basis) for each channel. 

From the eigenvalues, we create effective masses 



aMfj|(t + a/2) = ln 



A(i,io) 



in) 



A(i-f a,to)(") 



(34) 



Plotting these versus t, we look for plateaus of at least 
three points, within the single-elimination jackknife er- 
rors, when searching for a good fitting range in t. We 
also look for stable eigenvectors over the same region (the 
so-called "fingerprint" of the state; these are discussed in 
more detail in the next subsection, in the context of the 
couplings) . 

Figures [5] - 2] show samples of meson effective-mass 
plots from the two finer quenched ensembles and the 
larger dynamical lattice. The horizontal lines show the 
-|- Qr(7(") values from correlated fits (as before, we 
construct the covariance matrix with a single-elimination 
jackknife procedure) to the corresponding time ranges. 
Figure [S] shows similar plots for the baryons. 

We try to be more conservative with our choice of 
to here (previously we used to /a = 1) and work with 
to /a = 3 where we can [5l|. The notable exceptions are 



2. Extrapolations and interpolations 

Having the lattice energies, we may now reference 
each of them with respect to, say, the ground state: 
E - Eis = M - Mis + 0{l/m^). Figures M~E show 
these mass differences (in units of ro) as functions of the 
light-quark mass on the two finer, quenched lattices and 
the larger, dynamical lattice. (The pion masses associ- 
ated with the light-quark valence masses were computed 
from li ght- light correlators using standard source tech- 
niques j39l . [43 . [5^.) The experimental mass differences 



.(*) 



B 



(*) 



(*) 

sl(2) 
^(*) 



B. 



(*) 



Ab-B(*\ 



and — B^^^(^ are plotted as the 



s results at that value; 



green crosses (x) using ro = 0.49 fm (note that, for the 
last difference, we use the average mass of the and 
-Bs*'' plotted at {rriud + 'm-s)/2, enabling an easier com- 
parison with the Maq — Mi 
see Sec. IV D] for further discussion of this) 

For the most part, the plots display rather slight de- 
pendences of M — Mis on the light-quark mass and we 
extrapolate (interpolate) to mud {ms) by fitting linear 
functions of M^. For a few of the chiral (strange) ex- 
trapolations (interpolations), we leave out some of the 
points. These include the lightest quark mass for the 
25 state on the /3 = 7.57 lattice (due to a questionable 
rise there; see Table I VIII) and the heaviest quark mass 
for the Aq on the /? = 7.57 and (3 = 8.15 lattices (due 
to rather large values for ; see Tables |V] and IVIIP . We 
point out, however, that the corresponding extrapolation 
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FIG. 2: Sample of meson efTective masses from the 16^ x 32, f3 = 7.90 quenched lattice {aniq = 0.08). From left to right are 
the 5", -D±, P_, and P+. The horizontal lines denote the aE^'^^ :taa^'^^ correlated fit values and the corresponding time ranges. 
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□ X, 
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o X, 
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FIG. 3: Sample of meson effective masses from the 20^ x 40, /3 = 8.15 quenched lattice {arriq = 0.064). From left to right are 
the 5", D±, and P+. The horizontal lines denote the aE^'^^ ztctfj'-'^^ correlated fit values and the corresponding time ranges. 
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amq=-0.02, jS-4.65 


□ A, 
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P_, 16^X32 
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□ Xp 
I 



am =-0.02, ^=4.65 



□ xj 



t/a 



t/a 
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FIG. 4: Sample of meson effective masses from the 16^ x 32, (3 = 4.65 dynamical lattice (am, = —0.02). From left to right are 
the S, D±, P_, and P+. The horizontal lines denote the aE^"^ iaa'"' correlated fit values and the corresponding time ranges. 
The 25 states are extracted from two-mass fits to the eigenvalues A^'^-' with to/a = 1. 



results change by an amount within the final error bars 
when these points are left out. 

The results of the light-quark chiral extrapolations of 
the M — Mis A/ — Mg(,)) values for the lower-lying 

meson excitations and baryons (25', IP, ID, Af,, 
are displayed in Table |Xl Also included are continuum 
extrapolations of the quenched results, using linear fits 
in (see further discussion of this below). Since some 
of the fits for the continuum extrapolations, or some of 



the points which are used in them (labeled with ?), are 
rather poor (e.g., large x^'s), we present alternative fits 
for these in an attempt to check for strong systematics. 
For example, each of the eigenvalue fits for the 1P+ states 
on the /3 — 7.57 lattice have x^/dof > 2; so we also try 
a constant fit in a to just the two finer quenched lattices. 
The result, however, is only 1.2(T higher and both the 
old and new continuum-fit values are < 1 with one 
degree of freedom. 
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FIG. 5: Sample of baryon effective masses from the /3 = 7.90 and {3 — 8.15 quenched lattices and the [3 = 4.65 dynamical lattice 
(from left to right, respectively). The horizontal lines denote the aE^"^ zbaa'-"' correlated fit values and the corresponding time 
ranges. The ground states are extracted from two- mass fits to the eigenvalues A'^'. 



Contrary to our expectations and the results on all of for the 25* on the coarse quenched lattice lands above 
the other lattices, the result of the chiral extrapolation that for the 1D±. For this reason we try a constant 

continuum fit to the two finer lattices and find a result 



TABLE V: Results from correlated mass fits on the /3 = 8.15 
quenched lattice. The f symbols indicate where two-mass fits 
are used. 



TABLE VL Results from correlated mass fits on the (3 = 7.90 
quenched lattice. The f symbols indicate where two-mass fits 
are used. 



state 


aniq 


to/a , 


t/a 


aE 


x'/dof 


state 


arriq 


to /a , t/a 


aE 


x'/dof 


IS 


0.016 


3 , 6- 


-15 


0.4824(15) 


2.85/8 


IS 


0.020 


3 , 6-11 


0.5548(24) 


0.56/4 




0.032 






0.4935(12) 


1.51/8 




0.040 




0.5688(17) 


0.65/4 




0.064 






0.51612(93) 


1.97/8 




0.080 




0.5968(13) 


0.91/4 




0.080 






0.52783(85) 


2.31/8 




0.100 




0.6108(12) 


0.90/4 


2S 


0.016 


3,4- 


-11 


0.8837(99) 


0.67/6 


2S 


0.020 


3 , 4-8 


1.065(17) 


0.15/3 




0.032 






0.8881(86) 


0.56/6 




0.040 




1.064(14) 


0.06/3 




0.064 






0.9000(72) 


0.50/6 




0.080 




1.081(11) 


0.22/3 




0.080 






0.9075(65) 


0.35/6 




0.100 




1.090(10) 


0.10/3 


IP- 


0.016 


3 , 6- 


-13 


0.7302(53) 


3.70/6 


IP- 


0.020 


3 , 5-9 


0.8738(77) 


0.98/3 




0.032 






0.7371(41) 


3.26/6 




0.040 




0.8830(57) 


2.88/3 




0.064 






0.7549(34) 


4.55/6 




0.080 




0.9009(43) 


4.83/3 




0.080 






0.7649(31) 


5.56/6 




0.100 




0.9113(39) 


5.54/3 




0.016 


3 , 6- 


-13 


0.7425(49) 


2.58/6 




0.020 


3 , 6-10 


0.875(12) 


5.16/3 




0.032 






0.7472(37) 


4.07/6 




0.040 




0.8737(80) 


5.30/3 




0.064 






0.7634(28) 


6.26/6 




0.080 




0.8919(56) 


4.34/3 




0.080 






0.7727(26) 


7.08/6 




0.100 




0.9035(50) 


3.69/3 


1D± 


0.016 


3 , 6- 


-10 


0.930(19) 


0.77/3 




0.020 


3 , 5-8 


1.148(20) 


0.24/2 




0.032 






0.937(15) 


0.62/3 




0.040 




1.149(15) 


0.26/2 




0.064 






0.953(11) 


0.63/3 




0.080 




1.160(11) 


0.31/2 




0.080 






0.961(10) 


0.66/3 




0.100 




1.167(10) 


0.52/2 


Aq 


0.016 


3,4- 


-12 


0.7329(79) 


4.16/5+ 


Aq 


0.020 


3 , 4-11 


0.848(13) 


8.0/4+ 




0.032 






0.7644(47) 


6.08/5+ 




0.040 




0.8787(92) 


7.2/4+ 




0.064 


3 , 5- 


-13 


0.8262(44) 


13.4/5+ 




0.080 




0.9557(57) 


3.0/4+ 




0.080 






0.8586(38) 


15.7/5+ 




0.100 




0.9952(48) 


2.1/4+ 




0.016 


3,4- 


-12 


0.8521(66) 


4.05/5+ 




0.020 


3 , 4-9 


0.954(24) 


2.6/2+ 




0.032 






0.8711(43) 


2.42/5+ 




0.040 




1.0063(98) 


2.3/2+ 




0.064 






0.9152(30) 


3.70/5+ 




0.080 




1.0667(52) 0.78/2+ 




0.080 






0.9390(27) 


4.54/5+ 




0.100 




1.0959(43) 0.59/2+ 
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2Aa higher; however, both contimmm extrapolations 
have decent values. 

For the baryons on the finest lattice, we unfortunately 
have the smallest basis (thus far only 2x2) and are there- 
fore left with a suspicion that we may not have well iso- 
lated the ground-state mass. We thus try constant fits 
in a to the two coarser lattices and the results end up 
« 2.6(T and « 1.9a lower for the Af, and Ej, , respec- 
tively, with the continuum-fit /dof value for the latter 
greatly improving (from 4.4 to 1.8). This consistency of 
the mass splittings at the two larger values of a, together 
with the formerly poor Ej^ extrapolation, suggests that 
we may indeed need to enlarge the baryon basis on the 
P = 8.15 lattice to get a more reliable (perhaps lower) 
value of M (.) — Af^c.) (and maybe — M^(,) also). 

We also check the baryon masses on the coarser 
quenched lattices with a similar two-operator basis (only 
smeared). We find that the results change, but only 
slightly (i.e., well within the errors). This suggests that 
the higher-state terms on the fine lattice may also be well 
isolated, even with the small basis; but since the smeared 
operators are not quite the same on all lattices, we cannot 
be sure of this. 



TABLE VII: Results from correlated mass fits on the /3 — 7.57 
quenched lattice. The f symbols indicate where two-mass fits 
are used. 



state aniq to /a , t/a clE x^Jdof 

IS 0.025 3 , 5-10 0.7073(18) 3.21/4 

0.050 0.7227(15) 2.77/4 

0.100 0.7551(12) 2.02/4 

0.125 0.7715(11) 1.62/4 

2S 0.025 3 , 4-8 1.482(31) 4.67/3 

0.050 1.467(30) 4.97/3 

0.100 1.461(29) 0.20/3 

0.125 1.468(30) 0.07/3 

IP- 0.025 3 , 5-9 1.107(11) 1.05/3 

0.050 1.1184(86) 0.65/3 

0.100 1.1462(66) 0.63/3 

0.125 1.1605(60) 0.70/3 

1P+ 0.025 3 , 4-8 1.1569(46) 6.18/3 

0.050 1.1638(39) 6.92/3 

0.100 1.1843(31) 6.63/3 

0.125 1.1958(29) 6.38/3 

1D± 0.025 3 , 4-8 1.462(17) 0.51/3 

0.050 1.472(14) 0.84/3 

0.100 1.498(13) 1.37/3 

0.125 1.508(12) 1.34/3 

Aq 0.025 3 , 4-9 1.105(23) 3.01/2T 

0.050 1.163(12) 5.34/2t 

0.100 1.2609(80) 8.16/2t 

0.125 1.3061(71) 8.65/2^ 

E^*' 0.025 3 , 4-9 1.3027(53) 1.51/2^ 

0.050 1.3337(39) 1.42/2^ 

0.100 1.3936(66) 0.25/2^ 

0.125 1.4250(63) 0.05/2t 



Table IXII displays the strange-quark interpolations of 
the M — Mis {-^ M — M^(.)) values for the lower-lying 
states. Quenched continuum results are also present, 
again including some alternative fits for a few question- 
able cases: A similar IP+ constant continuum fit to the 
two smaller a values is displayed, resulting in a marked 
improvement in the fit (x^ = 5.4 — > 2.5), but an increase 
in the splitting of only « 1.6a. Although the coarse 
static-strange 25" falls below the 1D±, we nevertheless 
show the same constant continuum fit as for the static- 
light case and get a number « 1.8f7 higher, but again, 
with reasonable x^ values for both types of fits. As with 
the S^*-*, the rj^*"* continuum extrapolation is poor 
(x^ — 7.9 !) and the two results at larger a agree quite 
well, resulting in a significantly lower fij^ —Bg splitting 
(by K, A.Ta) with a constant fit to just these two. 

Figures [9] and [10] show the static-light and static- 
strange mass splittings as functions of the squared lattice 
spacing. The results at — > correspond to the lin- 
ear extrapolations using all three (a/ro)^ values. In Fig. 
[11] the baryon- meson mass differences are plotted versus 
(a/ro)^. Shown here are both the linear-a^ extrapola- 



TABLE VIII: Results from correlated mass fits on the /3 = 
4.65 dynamical lattice. The f symbols indicate where two- 
mass fits are used. 



state 


aniq 


to/a , t/a 


aE 




15 


-0.07 


3 , 5-12 


0.7328(30) 


7.81/6 




-0.06 




0.7420(26) 


6.57/6 




-0.05 




0.7508(24) 


5.79/6 




-0.02 




0.7772(20) 


3.47/6 


25 


-0.07 


1 , 2-9 


1.261(47) 


4.17/41" 




-0.06 




1.271(42) 


3.72/4t 




-0.05 




1.281(38) 


3.19/4t 




-0.02 




1.301(30) 


1.84/41 


1P_ 


-0.07 


3 , 4-9 


1.087(13) 


1.98/4 




-0.06 




1.087(12) 


1.33/4 




-0.05 




1.091(11) 


1.00/4 




-0.02 




1.1096(88) 


0.53/4 




-0.07 


3 , 4-8 


1.081(12) 


1.57/3 




-0.06 




1.083(11) 


1.80/3 




-0.05 




1.0870(95) 


1.99/3 




-0.02 




1.1065(78) 


2.61/3 




-0.07 


3 , 4-7 


1.430(56) 


1.14/2 




-0.06 




1.439(49) 


1.35/2 




-0.05 




1.445(45) 


1.53/2 




-0.02 




1.457(37) 


1.77/2 


Aq 


-0.07 


3 , 4-10 


1.054(39) 


O.34/3I 




-0.06 




1.078(30) 


0.25/3+ 




-0.05 




1.103(23) 


O.29/3I 




-0.02 




1.181(14) 


I.4O/3I 




-0.07 


3 , 4-9 


1.191(32) 


1.24/21 




-0.06 




1.205(27) 


1. 31/21 




-0.05 




1.221(22) 


I.I9/2I 




-0.02 




1.273(15) 


O.88/2I 
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tions and the constant fits to the two coarser lattices, 
the differences between the two sets being clear. 

As mentioned above, the baryon-meson mass differ- 
ences on the fine lattice are significantly larger than 
on the coarser lattices, and by so much that the 



continuum extrapolations for our Af„ 



M 



and 



M-g(,) values appear to be problematic. If we 



believe there are still some significant higher-order cor- 
rections in the fine baryon masses, we might try mitigat- 
ing some of these effects by considering differences of just 
baryon masses. Accordingly, in Table IXlTl we present re- 
sults of chiral (strange) extrapolations (interpolations) of 
M — values, along with quenched continuum extrap- 
olations. An associated plot of M — versus (a/rp)^ 
appears in Fig. 1121 These extrapolations appear much 
improved and suggest that we have to work harder to 
get reliable baryon-meson mass differences on the fine 
lattice. 

In Tables IXIIII and fXIVI we present results of chiral ex- 
trapolations and strange interpolations, respectively, of 
M — Mis rnass differences involving higher-lying states: 
35", 2P, 2D, A'ij. (Also included are quenched contin- 



TABLE IX: Results from correlated mass fits on the /3 = 5.20 
dynamical lattice. The f symbols indicate where two-mass fits 
are used. 



state 


aniq to /a 


, t/a 


aE 




IS' 


0.02 


3 , 


6-12 


0.6223(57) 


6.18/5 




0.04 






0.6390(45) 


4.39/5 




0.08 






0.6729(35) 


4.30/5 




0.10 






0.6904(32) 


4.90/5 


25' 


0.02 


3 , 


4-9 


1.024(82) 


0.96/2+ 




0.04 






1.038(66) 


0.56/2+ 




0.08 






1.060(64) 


0.28/2+ 




0.10 






1.074(64) 


0.24/2+ 


IP- 


0.02 


3 , 


6-10 


0.852(42) 


2.42/3 




0.04 






0.855(33) 


2.63/3 




0.08 






0.885(25) 


3.50/3 




0.10 






0.901(23) 


3.61/3 


1P+ 


0.02 


3 , 


6-10 


0.852(29) 


2.52/3 




0.04 






0.866(22) 


2.62/3 




0.08 






0.897(17) 


1.35/3 




0.10 






0.913(15) 


0.93/3 


1D± 


0.02 


3 , 


4-7 


1.114(62) 


0.26/2 




0.04 






1.125(53) 


0.20/2 




0.08 






1.157(44) 


0.37/2 




0.10 






1.173(41) 


0.46/2 


Aq 


0.02 


3 , 


4-9 


0.783(103) 0.21/2+ 




0.04 






0.893(50) 


0.61/2+ 




0.08 






1.024(27) 


1.23/2+ 




0.10 






1.076(22) 


1.21/2+ 




0.02 


3 , 


4-9 


1.001(32) 


1.05/2+ 




0.04 






1.027(23) 


1.39/2+ 




0.08 






1.101(16) 


2.80/2+ 




0.10 






1.139(15) 


3.15/2+ 




FIG. 6: Meson and baryon masses, relative to the meson 
ground state Mis = M{Bq*^), as a function of (oc iriq) 
on the 16'^ x 32, f3 = 7.90 quenched lattice. All masses result 
from fits to the eigenvalues of the complete 4x4 basis. Circles 
represent the S states; squares, the P_ states; diamonds, the 
P+ states; bursts, the 1D±; plus signs, the baryons Aq and 
E^'; and crosses, the experimental results [J 0, S Q S 



(usmg ro 



0.49 fm). 



uum extrapolations for the 2P states.) We consider these 
masses to be somewhat more dubious since they arise 
from fits over rather short t ranges (see Figs. [5] -[5]). 



B. Decay constant ratios 

In this subsection we turn to the analysis of the state 
couplings to the local operators, which provide us with 
the decay constants. 

Using the eigenvectors from the solutions of the gen- 
eralized eigenvalue problem, Eq. we create the ap- 
propriate ratios of projected correlator matrices, i?(i) 
(see Eq. P^ ). for the S-wa.ve states and the case where 
the open index i refers to the local operator. From these 
altered correlators, we find very similar effective-mass be- 
haviors as for the original eigenvalues and we perform fits 
to the i?(i)j-"-* over the same t ranges as before. From the 

amplitudes to these fits we achieve the quantity v["'^*v'l"^ . 
Assuming real overlaps, we arrive at the desired couplings 
(see Eq. (|15p ) and we can then form the appropriate ra- 
tios for, say, fs'/fB, (see Eq. ([11])) or /bJIb (see Eq. 

Figure [13] shows our resulting values for the static-light 
JbsI Ib versus the lattice spacing, a/r^. Both quenched 
and dynamical points are plotted, as well as results from 
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FIG. 7: Meson and baryon masses, relative to the meson 
ground state Mis = M{Bq*-'), as a function of on the 
20^ X 40, /3 = 8.15 quenched lattice. All meson masses result 
from fits to the eigenvalues of the complete 4x4 basis. The 
baryons masses are from the smaller 2x2 basis (smeared 
only) . The symbols have the same meaning as in Fig. [6] 
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FIG. 8: Meson and baryon masses, relative to the meson 
ground state Mis = M{Bq*^), as a function of M^ on the 
16^ X 32, f3 = 4.65 dynamical lattice. All meson masses result 
from fits to the eigenvalues of a 3 x 3 basis. The baryons 
masses are from the complete 4x4 basis. The symbols have 
the same meaning as in Fig. [6] 



the /3 = 7.90 lattice where we use the "naive" (or non- 
Hyp-smeared) links for the static quark (the green circle) 
and a smaller 3x3 basis (the magenta diamond) . One can 
easily see that each of these additional points agrees with 
the original result using the Hyp-smeared static quark 
and the 4x4 basis. 

As we are using improved gauge and fermion actions, 



TABLE X: Static-light mass differences and their quenched 
continuum values (linear fits in a^). Values are given using 
ro = 0.49(1) fm [4J|. The ** symbols indicate constant fits in 
a using only the two finer quenched lattices (thereby avoiding 
some possibly problematic fits on the f3 — 7.57 lattice). The 
f symbols indicate constant fits using only the two coarser 
quenched lattices (these are the only ones which use the full 
4x4 basis for the baryons). 



operator states /3 



M ~ M 



(MeV) 



25 



7.57 800(29)(?) 

7.90 695(24) 

8.15 683(18) 

oo 612(31)(13) (x^ = 0.67) 

oo(?) 687(14)(14) (x^ =0.14) 

4.65 674(66)(14) 

5.20 712(139)(15) 



1P_ 


Bli 7.57 


408(12) 




7.90 


438(11) 




8.15 


420(9) 
435(15)(9) (x' = 2.81) 




CO 




4.65 


454(19) (9) 




5.20 


398(78) (8) 



B 



FT 

1(2) 



7.57 465(6) (?) 

7.90 433(16) 

8.15 441(9) 

oo 423(13)(9) (x^ = 0.81) 

oo(?) 439(8)(9) (x^ = 0.20) ** 

4.65 446(17) (9) 

5.20 402(54) (8) 



ID± 


r(.',*) 
-°2(3) 


7.57 


771(18) 




7.90 


811(27) 






8.15 


759(33) 
794(42)(16) (x^ = 1.80) 






oo 






4.65 


896(78)(18) 






5.20 


858(115)(18) 


Aq 


Ab 


7.57 
7.90 


340(25) 
361(17) 






8.15 


389(12)(?) 
415(23)(8) (x^ = 0.41) 
355(14)(7) (x' = 0.49) f 






oo 






c»(?) 






4.65 


358(55) (7) 






5.20 


260(142)(5) 






7.57 
7.90 


574(4) 
545(21) 






8.15 


601(11)(?) 
604(16)(12) (x^ = 4.41) 
573(4)(12) (x' = 1.79) f 






oo 






c»(?) 






4.65 


555(47)(11) 






5.20 


611(58)(12) 
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we were originally expecting to use linear extrapolations 
in (J? for these decay constants. However, since we do not 
improve the local static-light currents, we check these 
results for significant 0{a) effects by extrapolating the 
quenched points (where we have three a values) using 
three different functional forms: constant, linear-a, and 
linear-a^ fits (the resulting extrapolations go from left 
to right, respectively, in the plot). Not surprisingly, the 
constant fit (which is dominated by the largest a value) is 
inconsistent with the other two extrapolations, showing 
that there are likely lattice-spacing effects in this ratio. 
The difference between the other two central values is 
almost as large as that between the constant and linear- 



TABLE XI: Static-strange mass differences and their 
quenched continuum values (linear fits in a^). Values are 
given using ro = 0.49(1) fm [Sj. The ** symbols indi- 
cate constant fits in a using only the two finer quenched lat- 
tices (thereby avoiding some possibly problematic fits on the 
/3 = 7.57 lattice). The f symbol indicates a constant fit us- 
ing only the two coarser quenched lattices (these are the only 
ones which use the full 4x4 basis for the baryons). 



operator 



IS 



states 



M-M^^.) (MeV) 



7.57 722(30) 

7.90 657(15) 

8.15 647(12) 

oo 604(26)(12) (x^ = 0.42) 

oo(?) 651(10)(12) (x^ =0.24) 

4.65 664(39)(13) 

5.20 691(125)(14) 



IP- B*o,i 7.57 396(7) 
7.90 413(6) 
8.15 403(6) 
oo 412(10)(8) (x^ = 2.29) 
4.65 421(12)(9) 
5.20 382(59) (8) 



B 



1(2) 



7.57 435(4) (?) 

7.90 402(8) 

8.15 417(5) 

oo 400(8) (8) (x^ = 5.39) 

oo(?) 413(4)(8) (x^ = 2.54) ** 

4.65 417(10)(9) 

5.20 396(40) (8) 



-°s2{3) 



FT 



7.57 750(13) 

7.90 764(16) 

8.15 735(20) 

oo 749(26)(15) (x^ = 

4.65 863(47)(18) 

5.20 850(93)(17) 



1.33) 



7.57 642(6) 

7.90 639(8) 

8.15 674(6)(?) 

oo 683(9)(14) (x^ = 7.87) 

c»(?) 641(5)(13) (x' = 0.13) t 

4.65 624(21)(13) 

5.20 694(42)(14) 
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FIG. 9: Quenched continuum extrapolations of excited B 
meson masses, relative to B^*\ 




FIG. 10: Quenched continuum extrapolations of excited Bs 
meson masses, relative to Bi*\ 



fits, but they do overlap, given the rather large error 
for the linear-a fit. Therefore, we do not believe that we 
are able to conclusively determine whether the leading- 
order effects in a are negligible. However, since the linear- 
a and linear-a^ results are consistent (even more so for 
the other ratios; see below) , we choose to report our final 
value(s) from the linear-a^ extrapolation(s) (see Table 
IXVl) . 

The result from the larger dynamical lattice is signifi- 
cantly higher than the quenched result at similar lattice 
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FIG. 11: Quenched continuum extrapolations of baryon 
masses, relative to B'*' (for A;, and E[,*') and (for fi^*'). 
The bursts denote constant extrapolations using only the two 
coarser lattices (the ones where we have results for baryons 
from the full 4x4 bases). 



spacing. The finer dynamical result is consistent with 
the quenched, but the volume is rather small (1.35 fm; 
perhaps this enhances /b more than /b^ , reducing the 
ratio). 

Turning now to the first-excited S'-waves, we form ra- 
tios of the form of Eq. (fTB|) . In Fig. [T3] we plot these 
quantities from the two finer quenched lattices and the 
larger dynamical lattice as functions of the light-quark 
mass. The strange-mass interpolations then give the de- 



TABLE XIL Chirally extrapolated and strange-interpolated 
baryon mass differences and their quenched continuum values 
(where possible). Values are given using ro — 0.49(1) fm [4^ . 

operator states /3 M - Ma^ (MeV) 



"A^ 



"Al 8l5" 



7.57 228(21) 
7.90 184(29) 
8.15 228(16) 

oo 200(27)(4) ix^ = 1.20) 
4.65 195(72)(4) 
5.20 353(156)(7) 



n 



ITT 



900(68) (18) 



7.57 391(25) 
7.90 364(21) 
8.15 365(12) 

oo 350(23)(7) (x^ = 0.12) 
4.65 352(55)(7) 
5.20 390(142) (8) 



FIG. 12: Quenched continuum extrapolations of the S^*^ — A, 
and nj,*' - Ab baryon mass differences. 



cay constant ratios fs'/fB^, which we plot versus the 
lattice spacing in Fig. [I5l 

The most striking feature here is that we obtain 
[f B'J Ibs) static > 1- Perhaps naively, we would have ex- 
pected this ratio to be < 1 as the numerator involves 
the local current coupling to the excited state, which one 
could easily imagine to be less than that for the ground 
state (e.g., in a nonrelativistic potential- model view, the 
normalized wavefunction at the origin is smaller for the 
first-excited state than that for the ground state). 



TABLE XIII: Static-light mass differences, and their 
quenched continuum values (where possible), involving states 
whose successful isolation is more dubious. Values are given 
using ro = 0.49(1) fm gj. 

operator states ^ M - Mg(.-, (MeV) 



3S 



B 



(*)" 



7.90 
8.15 



1198(91)(24) 
1232(65)(25) 



2P_ Bo,i 7.57 1105(57) 
7.90 1093(72) 
8.15 986(49) 
oo 941(78)(19) (x^ = 0.78) 
5.20 1195(67)(24) 



2P+ 



B 



FT" 

1(2) 



7.57 1062(35) 
7.90 933(83) 
8.15 906(53) 

oo 811(79)(17) (x^ = 0.05) 
4.65 1047(72) (21) 



2D± Bl 



"A^ 



8.15 



1401 (64) (29) 



A 



8.15 



1289(68) (26) 
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FIG. 13: Ratio of meson decay constants {f Bs / fs) static as 
a function of lattice spacing. The three quenched, continuum 
extrapolations of the 4 x 4-basis results correspond to (from 
left to right): a constant fit, a fit linear in a, a fit linear in a^. 
The green circle at a/ro ~ 0.3 (shifted to the right for clarity) 
corresponds to the adjacent quenched result, but with naive 
links for the static quark. 



We therefore try to test the robustness of this rather 
strange resuh by looking for some systematic effects: We 
first tried reducing the size of the basis (to 3 x 3 for 
quenched) to see if the ratio is significantly affected by 



TABLE XIV: Static-strange mass differences, and their 
quenched continuum values (where possible), involving states 
whose successful isolation is more dubious. Values are given 
using ro = 0.49(1) fm [IJ. 



operator states 


P 


M-M^(.-, (MeV) 


3S Bi*^" 


7.90 


1134(57)(23) 




8.15 


1181(39)(24) 



2P_ B^o.i 7.57 1016(37) 
7.90 959(47) 
8.15 937(36) 
oo 891(56)(18) (x^ = 0.001) 
4.65 851(89)(17) 
5.20 1126(64)(23) 



2P+ 



B 



FT" 

1(2) 



2D± B 



WW 

s2(3) 



7.57 1011(26) 
7.90 910(53) 
8.15 886(33) 

oo 811(50)(17) (x^ = 0.05) 
4.65 949(41)(19) 



the fact that we do not have a truly orthogonal basis (the 
different overlaps can thus "trade" their relative contri- 
butions without sacrificing their usefulness in isolating a 
given state [111)- What we find is that, whenever we get 
the same energy for the first-excited state (from fits to 

(2) 

R{t)l ) from the 3- and 4-dimcnsional bases, we achieve 
the same results for the v'l'^^ and v^/''^ values, and hence 
Ib'/Jbs: too. The smaller basis, however, can some- 
times lead to higher results for v,- ' (due to higher-order 
corrections) and we find up to ^ 15% higher results for 
fpg/fps on the l3 = 7.90 lattice. These come together 
with larger values for E^^^^ , so we know in such cases that 
the larger basis better isolates the first-excited state, giv- 
ing more reliable (and lower) results for this decay con- 
stant ratio. It therefore appears that, with the 4x4 
basis, we have better isolated the 2S state for a determi- 
nation of fpg. However, without an even larger basis, we 
cannot strictly rule out that there are still higher-state 
corrections enhancing this quantity. 

Another test we apply is to attempt the extraction 
of the same numbers via a different method: direct fit- 
ting of the local-local correlators. Since we want an ac- 
curate amplitude for the first-excited state from this, 
it is necessary to perform such a fit to a three-state 
ansatz (otherwise, higher-state contamination will enter 
strongly at the level we are after). This proves rather 
difficult with our data and we are only able to get said 
fits to converge on our finest quenched lattice {(3 = 8.15). 
The results, however, are consistent with those from the 
R{t)'f^\ For example, at aniq = 0.064 we arrive at the 

result fpg/fps — 1.416(143) from a three-state fit to 
the local-local correlator over the range of t/a = 7 — 16 
ix^/dof = 8.2/4), with energies E'^^^ and E^"^^ consistent 
with those found from the variational analysis. Although 
this is not a very good fit and the window of meaningful 
fits is small (starts at t/a — 6 and 8 result in a much 
larger "x^/dof and degenerate i?^^-' and E'^^'^ values, re- 
spectively), this is only k, 1.1a above the correspond- 
ing result from jackknifed fits to the ratios R(t)f^ and 



1.263(20). 



TABLE XV: Static-light decay constant ratios and the 
quenched continuum values (linear-a^ fits). Values are given 
using ro = 0.49(1) fm fill to set the physical rris point (see 

Sec. Him}. 

P ifSs/ fB)static {f B'J fSs) static { f B' / f B') static 



8.15 



1301(44)(27) 



7.57 
7.90 
8.15 

oo 

4.65 
5.20 



1.042(7) 
1.058(21) 
1.075(26) 
1.087(31) 
(x' = 0.07) 
1.108(29) 
1.089(41) 



1.309(75) 
1.237(15) 
1.259(25) 
1.240(58) 
(x' = 1-34) 
1.356(142) 
1.453(168) 



0.976(142) 

0.996(59) 

0.972(61) 

0.972(123) 

(x' = 0.07) 

1.089(259) 

1.026(128) 
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FIG. 14: Sample of strange- mass (linear) interpolations of the quantity f^g/fp's on the (3 — 7.90 and (3 = 8.15 quenched 
lattices and the /3 — 4.65 dynamical lattice (from left to right, respectively). 
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FIG. 15: Ratio of meson decay constants (/s/ / fB^) static as 
a function of lattice spacing. The symbols have the same 
meanings as in Fig. 1131 



We also switch to unsmeared links for the static quark 
and check the results again. A sample of such a ratio 
on the (3 = 7.90 lattice appears as the green circle at 
a/ro w 0.3 (shifted to the right for clarity) in Fig. [151 
This value is obviously consistent with the one from the 
Hyp-smeared static quark. 

Further discussion of these tests, as well as that of the 
other decay constant ratios we determine, appears in Sec. 



FIG. 16: Three- to two-point correlator ratios 7?.'^'^' (above) 
and roT?.'-^' (below) for the first-excited S'-wave; see Eqs. (|25|) 
and (|27l) . respectively. The plateaus in t and t' show that 
the exponential terms from the first-excited static-light state 
have been successfully isolated in the projected three-point 
functions and canceled by the projected two-point functions 
in the denominator. Displayed is the case of i — j for the local 
operator in the basis, but the results are the same (well within 
the errors) for other combinations of i and j in TZ^'^'^\ These 
results are from 80 of the /3 — 7.90 quenched configurations 
at aniq — 0.08, with only one of the two possible boundary 
crossings of the static-light correlators. 



C. Kinetic corrections 



Results for the ratio /b' / Jb' are also presented in Ta- 
ble |XVl All of these values are consistent with 1. 



Applying the techniques presented in Sec. Ill Cl we cre- 
ate ratios of eigenstate projections of three- and two- 
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FIG. 17: Comparison of the "traditional" {T{t,t')22/C{t)22) 
and the variational- method (7?.'^') ratios used for extracting 
the kinetic correction to the ground state. The results agree 
a.t t — t' = 6a, consistent with the starting value of our chosen 
fit range for the ground-state eigenvalue (see Table|Vl}. These 
results are from 80 of the (3 = 7.90 quenched configurations 
at aniq = 0.08, with only one of the two possible boundary 
crossings of the static-light correlators. 



point correlator matrices, where the eigenvectors from 
the two-point variational method are used (see Eqs. (PT|) . 
(P^ . and (12Z1)). These we check for independence 

from the t and t' variables and the corresponding, bare 
kinetic corrections ef"'" ) may then be found. 

Figure [TBI shows our first results for the ratios TZ^f^^ 
and TqIZ^^'^ (versus t and f ) from 80 of the quenched 
f3 = 7.90 configurations at anig = 0.08. The first ratio 

should approach the quantity e''^'^V(''^i^''*^j-^'); the sec- 
ond rQS^'^'^\ Each of these shows plateaus (within the 
errors) in t and t' , indicating the proper canceUation of 
the exponentials associated with the first-excited state. 
We then perform correlated fits to obtain the bare kinetic 
corrections to the 25* static-light state. 

One should be careful in applying this method by look- 
ing for plateaus in the time range suggested by the two- 
point variational problem. For example, Fig. [T7] shows 
a comparison of the variational-method ratio T?.*-^-* with 
the usual way of extracting the correction to the ground 
state: T{t,t')22/C{t)22, where the subscript 2 refers to 
the smeared operator. Since we fit the two-point eigen- 
values for £''^^^ starting at t/a = 6, we use t'/a = 6 for 
the time of the Q&Q interaction. One can clearly see 
in the plot that the two methods only agree when t — t' 
also becomes 6a, even though one might be tempted to 



FIG. 18: Unsubtracted, bare kinetic-energy corrections to 
the ground and first-excited S'-waves (80 configurations, one 
boundary crossing) and the ground P+ (50 configurations, one 
boundary) from the (3 = 7.90 lattice, plotted versus the light 
quark mass. 



choose a seemingly earlier plateau for TZ^^\ Also, at any 
given time slice, there is no apparent statistical advan- 
tage to the variational method for the ground state; in 
fact, the "traditional" approach reaches an earlier, con- 
sistent plateau and for this reason we use this method to 
obtain the e'-'^^^) values. The results for n = 2, 3, however, 
involve time ranges consistent with what would be cho- 
sen for the energy and coupling determinations (starting 
at ^0 + la)- 

Fits to the relevant ratios for the 35* and \Pj^ states are 
also made and we plot some of the results versus the light- 
quark mass in Fig. 1181 Linear chiral extrapolations and 
strange interpolations are then possible and differences 
of the results are taken to cancel the Eg value (due to 
divergent QD^Q operator; see Eq. p^ ) and arrive at 
more meaningful shifts to the static-light mass splittings: 
see Table IXVTl Of course, these are still bare quantities 
in that the QD^Q operator, although tadpole-improved, 
requires a renormalization factor Z. 

Assuming Z k, 1 for the moment, we can form some 
rough estimates of the energy shifts to the static-light 
mass splittings. Working with a value of « 4.2 GeV 
we find that the 25* — IS* mass differences should 
increase by « 40(13) and « 52(8) MeV at = to„^ 
and m^, respectively. For 35" — IS, these work out to 
w 120(80) and « 120(50) MeV, while for 1P+ - IS* we 
find sa 33(30) and w 38(14) MeV. A systematic dis- 
crepancy in the n = 1 values becomes obvious when 



we calculate the 0(1/771^) adjustment to 



(1) 



E. 



(1) 
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in going from the static limit to the B-meson system: 
(ei^'^^ - ^2^6 « 2.4(1.2) MeV. Using experimen- 
tal, spin-averaged, heavy-hght mass differences to set the 
strange-quark mass in Sec. Illl Dl we found 9.9 MeV for 
this shift. We need to boost our statistics here (and luck- 
ily, we can and will do so) and run at different lattice 
spacings to see if such a systematic difference persists. 
For the 0(1/771^) mixing of the IS* and 2S static- 



strange states we find e 



(2,1) 



0.26(5) GeV2. If we have 



the sign correct and we assume Eg « and Z « 1, this 
would contribute to an « 6% decrease in fs^ and an 
4% increase in /s^, increasing Jb'JJbs by ~ 10% 
and exacerbating the "dilemma" we face in expecting 
Ib' I < 1- However, without an actual value for £q 
and without the mixings with states of n' > 2, we cannot 
accurately obtain the overall shifts in the couplings <5w,|"'* 
(see Eq. (gni)- 

Simply put, we have much more work to do here in or- 
der to go beyond these rough estimations of the kinetic 
corrections to the static-light states. We are encouraged, 
however, by the ability of the variational method to iso- 
late the excited-state contributions in these three-point 
function calculations. We point out again (see Sec. Ill C|) . 
though, that this method could exhibit instabilities. In 
the limit of large statistics or with a poor basis of op- 
erators, one will have to add the interaction (or current 
insertion), along with the subsequent transfer to the sink 
(and, separately, the source), to the basis since it is an 
additional, linearly independent operator. 



TABLE XVI: Preliminary differences of bare heavy-quark 
kinetic energy corrections on the /3 = 7.90 (a = 0.145(3) fm) 
quenched lattice. 



difference 



GeV^ 



,(2,2) 



e[Y> 0.02(1) (?) 
^(2,2) 0.12(7) 



e 

4^.1 0.06(12) 

^(5,2) 



^(2,2) 
^{3,3) _ 
^{3,3) 

g(l.l) 

,(1.1) _ 



-ud 
-(1,1) 

- s 

-(1,1) 

"ud 

-(1,1) 

- s 

-(1,1) 



0.34(11) 
0.44(7) 
1.0(7) 
1.0(4) 
0.28(25) 
0.32(12) 



V. DISCUSSION 

We now discuss our results separately for the S-, P-, 
and £>- waves, and the baryons. In Table lXVlll we provide 
a summary listing of the lower-lying static-light meson 
and baryon mass splittings, in which we include all lat- 
tice spacings for the quenched continuum extrapolations 
(see Sec. IIV A2l for details and caveats); all experimental 



values [l|, 0, H, H, [lOl for spin-averaged singly b hadrons 
are also listed. The decay constant ratios and preliminary 
kinetic corrections have already been well summarized in 
Tables KV] and Evil respectively. 



A. S- waves 

The ground-state S'-wave states are well determined on 
all ensembles. The lattice energies (i?*-^-*) we extract for 
these are very precise and provide good points from which 
the energies of all other states can be referenced. (A fu- 
ture implementation of the procedures of Ref. [36| should 
allow us to obtain values for Eq, Z, and Eg, and hence 
TO J and absolute hadron masses.) The differences among 
the ground-state energies also act as good measures of 



(1) 



the relative light-quark masses and we use E', 
set the strange quark mass (see Sec. Illl D[) . 

The local couplings to the ground states are also quite 
precise and we obtain values for the decay constant ra- 
tio {Jb J Ib) static with statistical errors of a few % (see 
Table |^| with O(IOO) configurations (< 2% error for 
the bare Jb^ on the /3 = 8.15 lattice). Quenched results 
for this ratio [5^ [H^l ar e typica lly lower than those from 
dynamical calculations [53, [53,113' by « 5 — 10%. Our 



TABLE XVII: Summary of our static-light (-strange) mass 
splittings (in MeV), approximating the spin-averaged mass 
differences of lower-lying &-hadron states. Above appear re- 
sults for which there are experimental values (f only the 
fii, has been observed experimentally); below, those without. 
Our values are given using ro = 0.49(1) fm [4J] and all lattice 
spacings for the quenched continuum extrapolations (see Sec. 
IIV A 21 and Tables IXl and IXTI for associated caveats). 



difference 


Nf = 


Nf = 2 


experiment 








a = 0.156(3) fm 


[1> 2, 3,4, 5, 50] 








M^,sea = 461(6)(9) 




r(*) 




423(13)(9) 


446(17) (9) 


423(4) 


n(*) 
«sl(2) 




400(8) (8) 


417(10)(9) 


436(1) 


Aft- 




415(23)(8) 


358(55) (7) 


306(2) 






604(16)(12) 


555(47)(11) 


512(4) 






466(17)(10) 


426(37) (9) 


476(5) 






683(9)(14) 


624(21)(13) 


765(16)'^ 




-At 


200(27)(4) 


195(72) (4) 


206(4) 




-Ab 


95(17)(2) 


111(37)(2) 


170(5) 




-Kb 


340(23)(7) 


342(55) (7) 


545(16)t 


B(*)' 


-bw 


612(31)(13) 


674(66)(14) 








604(26)(12) 


664(39)(13) 








435(15)(9) 


454(19)(9) 




^sO,l 




412(10)(8) 


421(12)(9) 




= ('.*) 
"6 


-Kb 


272(23)(6) 


269(55) (5) 




= ('.*) 
"6 




173(20)(4) 


158(50) (3) 
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quenched value (static- light, — > continuum limit: 
fsj Ib = 1.087(31))_is consistent with that in ^] and 
w 2(7 below that in [53] ■ The result from our larger dy- 
namical lattice is « 2(j higher than the quenched result 
at similar lattice spacing. The finer dynamical lattice 
gives a value consistent with the quenched, but repre- 
sents a rather small physical volume (1.35 fm); it may 
be that the local coupling is therefore enhanced more for 
the lighter quarks, reducing the ratio fsJIs on such 
a small lattice. Results on fine dynamical configurations 
at larger volumes and lighter sea-quark masses are neces- 
sary to properly handle the systematics described above. 
A thorough treatment of the necessary renormalization 
factors [53| would also permit us to go beyond decay con- 
stant ratios. 

Without explicitly removing the divergences associated 
with the static quark (£'0, ^o)' absolute masses and ki- 
netic corrections elude us for the moment. Up to a renor- 
malization factor, differences of the kinetic corrections 
determine the ©(l/rrij,) shifts to the static-light mass 
splittings. Our first results for ei^'^^ — e^^j^^ underesti- 
mate the expected ©(l/m^) correction to the i?!*'' — B^*^ 
mass difference (« 10 MeV) and we need to do further 
nmning of the three-point functions T{t, t')ij to improve 
our statistics and go to finer lattice spacing. Our pre- 
liminary results for the e^^'^) values are not yet precise 
enough for us to say anything about the corrections to 
the ratio 

With our choice of operator bases, the variational 
method [H, [13, HE HE 113, Ell works quite well for find- 
ing excited static-light 5'-wave states. In particular, the 
first-excited S'-wave (25) energies are well determined. 

We find quenched, continuum values for the Bg — Bl' 
and - mass splittings of around 600 MeV, the 
a K, 0.15 fm dynamical results being about 10% higher, 
yet consistent with the quenched values at similar lattice 
spacing. The few quenched energies which we find for 
the 35 state lead to estimates of the Bj^J^ — B^^^ mass 
differences of « 1.1 - 1.2 GeV (see Tables Km] and |XlV|. 

Hadronic decays of the B^*'> (and i?!*'' ) should pro- 
ceed through a two-pion L = emission, B''*^7nr, as the 
« 600 MeV puts these states well above this threshold. 
However, on our ensembles the sea quarks are too heavy 
(lightest Mtt, sea ~ 460 MeV) for this channel to open. 
The same is true for the two-step decay through the 
broad Bq ^. (Single pion decays directly to the ground 
state are possible, but partially suppressed since these 
would be L = 1.) The relative masses of our excited 5- 
wave states therefore suffer systematic shifts due to the 
suppression of these decay modes. 

The discussion of strong decays begs the question of 
how they should be handled in possible future studies 
and we therefore pause briefly here to comment on the 
subject. As the dynamical quark masses are reduced, 
multiple-hadron states with the same symmetries as the 
intended hadron may become energetically favored: e.g., 
in a lattice search for the B' at near-physical pion masses 



and large volumes, the Btttt, L = and Bit, L = 1 scat- 
tering states should more readily dominate the 5-wave 
correlators than the B' at large times. In (partially) 
quenched simulations it is known that the variational 
method can be used to isolate unphysical "ghost" states 
[IE] , and hence there is reason to believe that it may also 
work to separate hadronic scattering states as well. The 
inclusion of interpolators which better project the mul- 
tiparticle states may also become necessary, along with 
a finite- volume study to identify the scattering states, 
or alternatively, an on-shell-mixing analysis to determine 
partial widths and mass shifts (see, e.g., Ref. [l^). How 
much of this laundry list we will be able to cover in fu- 
ture simulations remains uncertain, but we shall certainly 
keep it in mind. 

Using the methods outlined in Ref. |26j] and Sec. IIIBl 
we are also able to find rather precise values for the lo- 
cal coupling to the first-excited state. Surprisingly, the 
resulting ratios, {fB'/fBs)static, are greater than 1 (see 
Table 1^1. Results from the hght-hght ^ and heavy- 
heavy |60|] sectors, as well as considering the heavy-light 
system in the context of nonrelativistic potential models, 
suggest that this should not be so. After some careful 
testing (see Sec. IIVB[) . we are left with a few possible 
explanations for this seemingly odd result: 

1. Higher-state corrections: We cannot completely rule 
out the possibility that there are still higher-state correc- 
tions in the eigenvectors. We check results from larger 
t values and find noisier, yet consistent values for this 
ratio. However, we note that some results from a smaller 
basis give larger excited-state couplings. Whenever this 
happens, we also find significantly larger values for the 
energies, so we know that the larger basis better isolates 
this state. But without an even larger basis, we do not 
know for sure whether the measured first-excited energies 
and couplings would be reduced further. 

2. The local operator has little effect in the basis: It 
could be that the local operator is not having a large 
influence on the ability of the method to isolate the 25 
state. If this were true, the associated eigenvector com- 
ponent would most likely have a value dictated by the 
other operators in the basis as their components did the 
work of separating the excited state. As far as we can tell, 
this does not seem to be happening. The decay constant 
ratios are quite similar on all ensembles, even though the 
smeared operators in the bases are not designed to be 
the same physical size. Results for the effective masses 
are similar whether we reduce the basis by removing the 
local operator or one of the others. Also, the few pos- 
sible three-state fits that we are able to perform on the 
fine, quenched, local-local correlators give similar num- 
bers for the ratio. However, such an effect as described 
above need only take place at about the 20 — 30% level 
(assuming an actual value of fs' / fn^ ^ 1) as we find 
[fsj fss) static ~ 1.2—1.3. More extensive testing of this 
method against direct fitting in bottomonia systems [2^ 
found no such systematic discrepancies. However, an- 
other group [27j . using a slightly different approach, did 
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see systematic enhancements of excited-state couplings 
in charmonia (e.g., /^'). 

3. The static approximation: It may be that the values 
for {f B'J fs,) static are correct, but are > 1 due to the 
use of the static approximation for the heavy quark. Our 
first step at relaxing this approximation (through £s ' 
at a w 0.15 fm), however, seems to increase the ratio 
further, but we reiterate here that without a value for £q 
and without the static-light mixings involving the second- 
and higher-excited states, we cannot accurately say. 

Of course, we are also left with the possibility that 
this effect is real for bottom-light systems. Indeed, some 
potential models which yield covariance display slightly 
higher decay constants for radially excited heavy-light 
states: in Ref. ^ f y/W / f^/B_^l.m is reported, 
whereas we find fB'^^/M^JfB,y/M^, = 1.31(6) in the 
quenched continuum hmit and 1.44(15) on the larger dy- 
namical lattice. 

Fortunately, for the corrections to mass differences we 
are in a position to say something more. The e(2,2)_^(i,i) 
differences thus far seem to indicate that the B^*^ — B^*) 
and Bl ' — Be mass splittings should be « 40 and ~ 
50 MeV higher, respectively, than the static-light values 
found on the quenched a « 0.15 fm lattice. 

B. P- waves 

The ground-state P+ states correspond to the spin- 
averages of the Bl (Bsi) and B2 (5^2 ) niesons recently 
observed at Fermilab 0, i. As can be seen in Table 
IXVIIl our results for the s| (2) ~ ^'*'' mass splitting are in 
reasonable agreement with the experimental value. How- 
ever, the true heavy-strange splitting, ^si(2) — Bs , hes 
above our determination. Note that the experimental 
heavy-light splitting has a lower value than the heavy- 
strange one, as expected [G^I- However, we see the wrong 
ordering here, due to the fact that the static results do 
not include the kinetic-energy correction for the heavy 
quark (using NRQCD up to such order [T^], the proper 
ordering was seen). We imagine that this will be resolved 
in the present approach when we have results for the ki- 
netic corrections from our full statistics (these should be 
larger for the heavy-strange mesons). What may then 
be different, however, is the overall, relative proximity 
of our and the experimental 1P+ — 15 splittings as our 
first results for the kinetic corrections (see Table IXVip 
suggest that these static-strange (and -light) mass differ- 
ences may increase by ~ 40 MeV (assuming mj, « 4.2 
GeV). 

As single-pion decays of the ^1(2) to the ground-state 
multiplet take place through L = 2 emission, these are 
relatively narrow resonances and the associated system- 
atic mass shifts should be small. The same is not the 
case for the Bq^^ (and i?*o_i)- L = 0, single-pion decays 
of these mesons are possible (see, e.g., Ref. [22|), making 
them rather broad and therefore difficult to see in exper- 



iments. Our 1P_ — IS* results (e.g., 1 ~ B^*'^) should 
thus only be viewed as "ballpark" figures since they likely 
include strong systematic effects due to the suppression 
of these modes. 

Calculation of the decay constants for the P-wave 
states is certainly possible [63|, and in the works here. 

Excited P-wave states are also observed, the mass 
sphttings with the 15 states lying in the 800- 900 MeV 
range (see Tables IXIIII and IXIV[) . The results here are 
somewhat less certain as the effective-mass plateaus are 
usually quite short. Oddly enough, the 2P_, states lie 
slightly above the 2P+ ones. More statistics and finer 
lattices (along with cross correlations with decay prod- 
ucts) are needed to better resolve these issues. 



C. D- waves 

We view the static-light D± meson as an approxima- 
tion of the spin average of the B* , B2, and P3 states. 
The reason for the presence of the B* in this average, 
rather than the B*, is the fact that the heavy-quark 
spin interactions which normally allow the mixing of the 
(mostly orbitally excited) 1D± configuration and the 15 
one are not present. Of course, the 1D± — 25 configu- 
ration mixing is therefore also absent; but these orbital 
and radial excitations are certainly closer in energy. The 
results appear to support this averaging scheme since the 
ID± energies are slightly above the level of the 25 on all 
lattices. 

All of the physical states associated with these op- 
erators are rather high-lying ones and should therefore 
experience a number of hadronic decays. Between our 
relatively heavy sea-quark masses and finite volumes (re- 
stricting the available momenta of decay products), we 
likely suppress these and we therefore have rough values 
for the mass differences here. 



D. Baryons 

With a four-operator basis on most ensembles and two- 
state fits to the largest eigenvalues of the variational 
method, we have a number of good results for the ground- 
state baryons. Differences between the resulting energies 
and those of the ground-state mesons are calculated, as 
well as those among the baryons (see Table IXVIIl) . 

There are obvious discrepancies between our quenched, 
continuum results for the Af, — P'-*^ and S^*'* — P^*-* 
mass splittings and the corresponding experimental val- 
ues, which are « 100 MeV (or « 5a) lower. Even consid- 
ering the possible systematic increase in the fine lattice 
results (due to only having a 2 x 2 basis there) and leav- 
ing them out, the baryon masses only decrease by « 60 
MeV for the Af, and « 30 MeV for the E^*^ (see Table 
|X|) . while the errors also decrease, thereby maintaining a 
largely significant difference with experiment. 
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The large- volume, dynamical results for these baryon- 
meson mass differences are consistent with experiment, 
but with much larger errors, making them also consistent 
with the quenched results at similar lattice spacing. 

The e[*^ - Ab mass differences (both quenched contin- 
uum and finite-a dynamical) agree with the experimental 

value. The fij, — Af, splitting also shows better agree- 
ment when comparing quenched continuum and dynam- 
ical results (but is below experiment; see below). This 
seems to support the notion that by taking differences 
between the baryon energies we are able to remove some 
of the possible higher-state corrections associated with 
the smaller basis on the fine quenched lattice (see Sec. 

IIV A2|) . Taking the s[*^ - A^ difference also partially 
removes the need for the heavy-quark relativistic correc- 
tions (see below) since the light valence content is the 
same for these two states. 

Given the presence of two lighter valence quarks in 
the baryons, the heavy quark in these states most likely 
exhibits a larger kinetic energy than in the mesons. Ac- 
cordingly, the kinetic corrections for our baryons may 
be quite large when compared to those for the ground- 
state mesons and the corresponding baryon-meson mass 
differences given here may become significantly greater. 
This could help to explain the much lower result we have 
for the nl*'^ - Bi*'^ mass splitting: « 80 MeV (« 9a) 

below the experimental fJj, — Bl for the quenched con- 
tinuum and K, 140 MeV (« 7a) lower for the dynamical 

a = 0.156(3) fm lattice. Also, since the il^,*'' has two 
strange quarks, while the Ab has two light ones, there 
is likely a sizeable amount of heavy-quark kinetic energy 
missing in our —A;, difference; indeed, our values are 
« 200 MeV too low. 

Hadronic decays and the associated mass shifts are also 
a more important consideration for the baryons and we 
have yet to consider any of these. 

Until now, we have not said much about the Sf, and 

' states. We do not yet explicitly create the appro- 
priate correlators for these as they require strange-light 
diquarks and we have only degenerate masses for ours. 
Nevertheless, we make the simplification of using our Aq 
results at = {m^^ + to approximate the Sf, 

and the corresponding Sg results for approximating the 

s{; '*\ These lead to the mass differences in Table IXVlII 
The only agreement with the experimental values comes 
from the quenched continuum limit of the Sf, — i?^*^ split- 
ting. The other values are all lower than experiment; the 
Sb — Ab mass difference being > 4ct away. This fits with 
our low — Ab observations due to the missing ex- 
tra heavy-quark kinetic energy for the baryons involving 
strange quarks. Given the possible over-simplification we 

use here for approximating the ^b and .^j, ' correlators, 
it becomes difficult to draw more meaning from these 
results. 

A few effective-mass plateaus for a possible, first- 



excited state are seen in the second Aq eigenvalues on 
the fine quenched (/3 = 8.15) lattice. The chiral extrapo- 
lation of these energies leads to an estimate of a A'^ — Ab 
mass difference of « 900 MeV (see Table |XII| . This rep- 
resents a large gap to the first excitation and, if real, 
would entail large systematic effects due to omitted ki- 
netic corrections and decay channels. We are, however, 
rather skeptical of this "state" as it arises in the case 
where we have the smallest basis (2 equally smeared op- 
erators with different diquark spin structures; see Table 
We see similar, yet somewhat less certain, plateaus 
on the /3 — 7.90 quenched lattice when we reduce the 
basis to the two-dimensional one similar to that above; 
these lead to a A'j, — Ab splitting estimate of « 970 MeV. 
However, a telling result is that the evidence of this state 
vanishes when we then enlarge the basis on this ensem- 
ble. Perhaps the same will happen with a larger basis on 
the fine quenched lattice. 



VI. CONCLUSIONS AND OUTLOOK 

In this work we presented static-light approximations 
of B- and B^-meson mass splittings and decay constant 
ratios in the quenched continuum limit and at two finite 
lattice spacings which include two flavors of improved 
sea quarks. Results were reported for a number of ex- 
citations: 25, 35*, IP, 2P, ID, 2D. We also included 
our first results for kinetic corrections to the IS, 2S, 3S, 
and 1P+ mesons, indicating the approximate shifts to 
the associated mass differences and decay constants. In 
addition, static-light-light baryon correlators and the cor- 
responding mass splittings were calculated. 

The strengths of the current calculation include: the 
use of improved gauge (Liischer-Weisz) and fermion (Chi- 
rally Improved) actions for the light degrees of freedom; 
a good separation of the static-light states, achieved via 
the variational method and a robust basis of hadron op- 
erators; and a significant boost in the statistics due to 
the domain-decomposition improvement (half- to- half) of 
the light-quark propagator estimation. 

Comparison of our results with experiment [l|, 0, S 0j 

fand with similar results from other lattice groups 
i, S [M [H Q [H, [13, M M M M is diflicuit 

for a number of reasons: mainly systematic effects inher- 
ent in our lattice formulation of HQET for the former, 
and those between the different formulations of HQET 
or QCD for the latter. However, we have pointed out 
the relevant literature above and have discussed our find- 
ings in detail. Overall we have shown: While our mass 
splittings from the ground state i?'*^ to the excited me- 
son states are quite reasonable, those up to the baryons 
are in general too large when near the light-quark chiral 
limit and too small for bottom-strange baryons. We dis- 
cuss possible sources for these discrepancies. Our decay 
constant for B'^ turns out to be somewhat larger than 
that for Bg. As this was initially unexpected, we per- 
form multiple tests of this result, but nevertheless find 
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it to persist. We offer a few scenarios whereby this may 
still be due to systematic effects, which we plan to test 
further. However, we also entertain the possibility of this 
being a real effect for heavy-light hadronic systems (see 
also, e.g., Ref. M)- 

Along the way, we have pointed out areas for further 
improvement. Higher statistics, finer lattices, larger vol- 
umes, lighter sea quarks, and more basis operators have 
all been mentioned and we shall endeavor to boost each 
of these. In the near future, we aim to remove the di- 
vergences associated with the static quark, make a de- 
termination of the bottom-quark mass, and calculate the 
kinetic corrections for all of the lower-lying hadrons on 
the full range of ensembles contained herein. 
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